Behavior of Phantom Scalar Fields near Black Holes 
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We present the accretion of a phantom scalar field into a black hole for various scalar field 
potentials in the full non-linear regime. Our results are based on the use of numerical methods and 
show that for all the cases studied the black hole's apparent horizon mass decreases. We explore 
a particular subset of the parameter space and from our results we conclude that this is a very 
efficient black hole shrinking process because the time scales of the area reduction of the horizon are 
short. We show that the radial equation of state of the scalar field depends strongly on the space 
and time, with the condition uj = p/p > —1, as opposed to a phantom fluid at cosmic scales that 
allows uj < — 1. 
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I. INTRODUCTION 

In cosmological models at present time the universe 
is assumed to have various ingredients of exotic nature 
like the dark energy candidates, among which we find the 
phantom scalar field pj. The reason is that supernovae 
la data allow an equation of state of the dark energy 
component with uj — p/p < —1, where p and p are the 
pressure and the energy density of the fluid, which is a 
condition that a phantom scalar field satisfies at cosmic 
scales. This is a reason to start up an exploration of the 
consequences such type of scalar field might have at local 
scales, for instance on black holes. On the other hand, 
the behavior of the area of the horizon when the black 
holes is accreting matter is a very important property of 
black hole physics, because the accretion of such exotic 
material may impose important restrictions on the mass 
of black nowadays black hole candidates. 

In order to study this process in the full non-linear 
regime we start up with a model coupling the phan- 
tom scalar field and gravity. The model assumes the 
Lagrangian density of the phantom scalar field is given 
by 



C 



(i) 



where R is the Ricci scalar of the space-time, g^y is the 
space-time metric, is a scalar field and V{cf>) is the 
potential of the field. The property defining a phan- 
tom scalar field is that the relative sign between the 
Ricci scalar and the kinetic term are the same. When 
the action constructed with such a Lagrangian density is 
varied with respect to the metric, the arising Einstein's 
equations are related to a stress energy-tensor that vi- 
olates the null energy condition, that is T^ v k^k 1 ' < 0, 
where fc M is a null vector. The immediate implication of 
this violation is the violation of the weak energy condi- 
tion too, which in turn implies that observers following 
time- like trajectories might measure negative energy den- 
sities. Although this property is at odds with the nowa- 



days physics observed in laboratories, there are cosmo- 
logical observations indicating the presence of such kind 
of matter In this paper we explore the implications of 
this at astrophysical scales. 

The fact that the scalar field violates the null energy 
condition motivates the study of possible unusual impli- 
cations in astrophysical scenarios, because the area in- 
creasing theorem does not apply in this case. For in- 
stance, using exact solutions corresponding to stationary 
accretion of a test phantom fluid, it was found that the 
mass of black holes decreases in a phantom energy dom- 
inated universe approaching the big rip Q. It has also 
been studied recently the behavior of a black hole appar- 
ent horizon (AH) in a FRW background, and the condi- 
tions under which a naked singularity can be formed due 
to the coincidence of the AH of the black hole and the 
cosmic horizon Q ; in such case the authors consider the 
effects of the back reaction of the scalar field on the space- 
time metric and indicate that cosmic censorship does not 
only forbid the existence of naked singularities but also 
the existence of a phantom field. Assuming that the null 
energy condition is satisfied the area of the horizon only 
increase. However, if this condition is violated, the area 
of the event horizon decreases and the black hole shrinks 
as shown recently in the full non-linear regime using nu- 
merical relativity in The present paper is a detailed 
follow up of previous one, in which we also explore the 
effects of the scalar field potential on the accretion rates 
and final state of the black hole. 

Two important items are presented in this paper: i) 
the relation p/p < —1 is not fulfilled at local scales by 
the scalar field (at least near to a black hole) although the 
null energy condition is not satisfied and ii) the accretion 
of such scalar field reduces the area of a black hole at 
similar rates for different types of potentials driving the 
scalar field. 

This paper is organized as follows. In section |H] wc 
describe the 3+1 decomposition of the space-time, the 
evolution system of equations driving the evolution of 
the geometry and the construction of the initial data. In 
section IIIII we present the results obtained. Finally in 



section HVl wc draw some conclusions and comments. A. Evolution 



II. THE SYSTEM OF EQUATIONS 

We formulate Einstein's field equations coupled to the 
scalar field in such a way that these can be integrated nu- 
merically. In this model, we will use spherical symmetry 
and use geometrized units for which the speed of light 
and Newton's constant are equal to one. 

In order to the study the dynamics of a spherically 
symmetric space-time, we write the general metric for 
the coordinate system (t, r, 8, #) in the following form 

ds 2 = -(a 2 - ^/3 r ^L)dt 2 + 2p r dtdr 
X 

+ - [g rr dr 2 + g 9 e{d9 2 + sin 2 9d(f> 2 )} , (2) 
X 

where x ac ts as a conformal factor relating this metric to 
a space-like flat metric, (3 r is the only non-zero compo- 
nent of the shift vector and a is the lapse function as in 

!• 

We solve the evolution Einstein's field equations using 
the Generalized BSSN evolution formulation of the 3+1 
decomposition of General Relativity described in [f| as 
opposed to previous successful analyzes refereed to the 
accretion of scalar field using Eddington-Finkelstein like 
coordinates under the ADM formulation Q. 



For the construction of the space-time we carry out a 
Cauchy-type evolution of initial data based on the 3+1 
decomposition of the space-time. Within such decom- 
position there are various formulations of the evolution 
equations with different hyperbolic properties. Among 
such formulations the most popular nowadays is the so 
called BSSN formulation Q which helped at solving the 
problem of the binary black hole collision system recently 
[e| Q using the punctures technique that allowed the ad- 
equate treatment of the black hole singularity [l(J. The 
BSSN formulation assumes that the conformal metric has 
determinant equal to one, but in spherical coordinates 
the flat metric the determinant is different to one. This 
issue is addressed relaxing this condition over the deter- 
minant, obtaining the Generalized BSSN equations (GB- 
SSN) [ll| , which is the formulation we use in our simula- 
tions. The GBSSN system of equations in a spherically 
symmetric space-time reduce to a set of six equations 
for the independent dynamical variables X i 9rr, gee, 

the 

non-zero trace- free part of the extrinsic curvature A rr , 
the trace of the extrinsic curvature K and the contracted 
Christoffel non-zero symbol (conformal connection func- 
tion) r r . The explicit expressions in the presence of mat- 
ter are: 



d t g rr = -2A rr a - ^ + P% r - ^EJk + 2 9r ^' - 2 -*^K 
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where primes denote derivatives with respect to r and dinates are Eulerian, whereas for v = 1 the coordinates 
the gauge parameter v is such that for v = the coor- are Lagrangian. 
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Wc have introduced the matter sources p, S r , M rr and 
S which represent the energy density, the momentum cur- 
rent density, the non-zero trace free part of the stress 
tensor and the trace of the stress tensor respectively. 
These quantities can be computed after projecting the 
stress-energy tensor onto the unit normal vector n a to 
the space-like hyper-surfaces, 



tests. In all our simulations we use Eulerian coordinates, 
that is, we set v = 1. 

The evolution of the scalar field is driven by the Klein- 
Gordon equation 



dV 



--d li [^g^d^]- — =0. (6) 



d(t> 



p = n a npT aP , 

S r = -lr a npT a P, 
M rr = lralrpT afi 

The constraints of the system are: the Hamiltonian 
constraint H, the momentum constraint Mi and the con- 
straint arising from the definition of the conformal con- 
nection functions G % . In spherical symmetry they are: 



H = - 



mi 
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(3b) 
(3c) 



In order to specify the gauge, we evolve the lapse ac- 
cording to the standard 1+log slicing condition 



d t a = f3 a d a a - 2aK, 



(4) 



or the related condition obtained by dropping the advec- 
tion term f3 a d a a. For the shift vector we implemented 
the recipe for the T-driver condition 



d t /3 a = -B a + (3 c d c /3 a , 

d t B a = d t T a + p c d c B a - [3 c d c T a - i]B a , 



(5) 



which helps to avoid the slice stretching near the horizon, 
which is known to kill the numerical evolution. In order 
to avoid instabilities near the origin (the puncture), we 
implement a sort of excision without excision [l2| using 
a factor function on the source of the evolution equations 
of the form (r/(l + r)) 4 from the coordinate origin out 
to r ex an appropriate fraction of the size of the apparent 
horizon radius such that r ex is smaller that the apparent 
horizon of the final black hole. Even though this function 
violates the constrains, the violations do not propagate 
out of the black hole horizon as shown by the convergence 



In order to reduce this equation to a first order system as 
the geometric counterpart, we define two new variables 
7r = dt(j> and ip = 9 r 0, and ([6]) reduces to the system 



dt-K = d t (3 r + (-7T + (3 r )d t \n 



ax 
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+ d r In 
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which together with the evolution equations for the ge- 
ometry and the completes the system of equations to be 
solved. 



B. Initial data 

In order to start up an evolution including the mat- 
ter terms it is necessary to solve the constraints a the 
initial time slice, and afterwards Bianchi identities guar- 
antee that such constraints are satisfied. We assume that 
the initial slice is time-symmetric, which implies that the 
shift vector j3 r and its time derivative are zero initially. 
In addition, time symmetry imply that all the compo- 
nents of the extrinsic curvature are zero initially. In this 
way, the momentum constraint is satisfied identically at 
initial time. On the other hand, we use the pre-collapsed 
condition for the lapse a = (1 + M/2r)~ 2 . 

In order to provide a local nature to the scalar field 
we assume that <f> has a Gaussian profile or is a train 
of Gaussians one next to the other, which implies initial 
data for the first order variable ip, and the time symme- 
try provides the condition tt = 0. With this information 
about the matter source we solve the Hamiltonian con- 
straint. Now, in order to achieve a smooth coordinate 
system we use the initial ansatz that the space-time met- 
ric has a form similar to that of a Schwarzschild black 
hole in isotropic coordinates. Thus we assume that the 
metric at initial time has the form 



9rr = 1, 



1 + ^ +U 



(8) 



where M is the mass of the apparent horizon and u = 
u(r, t = 0) is the function to be determined through the 
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solution of the Hamiltonian constraint. The Hamiltonian 
constraint is thus reduced to the equation 



2d r u 



v{4>) 



1 M 
1+ 27 



I 



M 
2r 



(9) 



which we solve using a fourth order Runge-kutta integra- 
tor. In the last equation ([9]), V(4>) is the self- interaction 
potential of the phantom scalar held. 

We finally rescaled the coordinate r in such way that 
X — > 1 when r — > oo. These initial data are used to start 
up a Cauchy type evolution using the GBSSN evolution 
equations (O, the gauge conditions dU-O and the scalar 
field equation ([7]). 



use this function in order to compare the mass of the 
apparent horizon and the mass of the space-time, so that 
we can monitor the consistency of our simulations. The 
Misner-Sharp mass function in our case is written as 



Mms 



R 



(d t R) 2 2p r {d t R){d r R) 



{d r Rf 



(12) 



where R = \Jgeejx- 

The Misner-Sharp mass reduces to the Arnowitt- 
Deser-Misner (ADM) at space-like infinity, 



M A 



DM 



lim M 



MS, 



(13) 



C. Implementation and diagnostics 

The numerical method used to approximate the con- 
straint and evolution equations is a fourth order finite dif- 
ferences approximation. We only perform the evolution 
on a finite domain with artificial boundaries at a finite 
value of r where we implement radiative-type boundary 
conditions. The integration in time uses a method of 
lines with a fourth order accurate Rungc-Kutta integra- 
tor. Throughout the evolution, we monitor the Hamilto- 
nian, Momentum and G r constraints, so that we check 
that they converge to zero as resolution is increased with 
fourth order. 

Apparent horizons. In order to track the radius, area 
and mass of the apparent horizon during the evolution in 
terms of the 3 + 1 variables, we search for the marginally 
trapped surfaces (MTS) through the condition 



so we can estimate the mass of the space-time at a fi- 
nite radius at any time during the evolution. On the 
other hand, the Misner-Sharp mass function reduces to 
the Bondi-Sachs mass at future null infinity. 

Event horizon. Aside of the apparent horizon which is 
a gauge dependent 2-surface, we track a bundle of radial 
outgoing null rays in order to approximately locate the 
event horizon during the evolution as the 3-surface from 
which outgoing null rays diverge when launched toward 
the future. For this we solve the geodesic equation for 
radial null rays on the fly during the evolution, and fine- 
tune their location initially so that both, outgoing null 
rays escaping toward future null infinity and outgoing 
null rays that reach the singularity converge to the same 
surface for as much time as possible during the evolution, 
after which they spread. Such surface is the approximate 
location of the event horizon. 



6 = V^ 4 + KijuW - K = 0, (10) 

where we take n l to be the outward pointing unit vec- 
tor normal to the horizon, Kij are the components of 
the extrinsic curvature and K its trace of a space-like 
hypcrsurface on which one calculates the MTSs. The ap- 
parent horizon is the outermost MTS. In our coordinates 
equation (fpQ|) for the metric ((2]) reads 

d r gee _ d r x _ [~g~^ ( K A rr \ 
9oo X V X U 2 9rr )-°- (H) 

In order to track the evolution of the apparent horizon we 
calculate O at every time step and locate the outermost 
zero of it at the coordinate radius tah- Then calculate 
the area of the corresponding 2-sphere Aah = 4nR\ H 
and its mass Mah = Rah I '2, where Rah = ^JWe ~fx * s 
the areal radius evaluated at tah- 

Misner-Sharp mass function This quantity is a space- 
dependent measure of the mass of the space-time. We 



III. RESULTS 

In Fig. [T]we show a particular case of an initial scalar 
field profile. We can see that the solution of the Hamilto- 
nian constraint implies already at initial time a negative 
energy density as expected from the violation of the null 
energy condition of this type of field. This also opens 
another possibility to be explored, that is, once we know 
that the scalar field contributes with a negative energy 
density it is also possible to construct initial configura- 
tions whose ADM mass is negative, that is, the contribu- 
tion of the scalar field energy density to the space-time is 
bigger in absolute value than that due to the black hole's 
horizon. 

On the other hand, we also show in Fig. Q]tlic relation 
lu = p/p where p = T r r , which clearly indicates that at 
local scales, the Lagrangian of the phantom scalar field 
does not behave at all as phantom matter since p/p > — 1, 
and in fact runs from cosmological constant uj = — 1 to 
stiff matter lu = 1. 
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FIG. 1. In the top panels we show the initial profile of <j> as function of the radial coordinate (in units of the initial mass of the 
apparent horizon) for the exponential potential in the case of positive ADM mass (a single Gaussian) and negative ADM mass 
(train of three Gaussians); these profiles were used to solve the Hamiltonian constraint at initial time. In the bottom panels we 
show the energy density of the scalar field p = T^n^n" labeled with the left axes, and the equation of state uj = p/p labeled 
with the right axes as function of the radial coordinate (using the same units as before) at initial time for two cases: i) the initial 
scalar field profile corresponds to a Gaussian (f>{r,t = 0) = 0.28e~' r_5 ' and exponential potential, and ii) the initial scalar field 

profile corresponds to the superposition of three gaussians <j>(r,t = 0) = 0.25e~' r ~ 5 ' + 0.125e~' r ~ 10 ' + 0.05e~' r ~ 15 ' , also for 
the exponential potential. At further values of time, the behavior of ui is the same, except that it appears shifted in the spatial 
direction. This behavior indicates that ui does not only depend highly on the space, but also on time. 



In our study we set up two classes of initial data cor- 
responding to positive ADM mass and negative ADM 
mass for all the potentials studied. In order to re- 
duce the parameter space of our study, for each poten- 
tial we fine-tune the initial scalar field profile so that 
the ADM mass has the same value for all the poten- 
tials explored. The potential we used are: a) quadratic: 
V((f>) = V cj) 2 , b) quartic: V{<j>) = V^ 2 + and c) 

exponential: V(cf>) = Vbe~*, where Vq and V\ are nega- 
tive constant parameters. We also explored the poten- 
tials V((f>) = Vb(cosh(/) — 1) and Vq = sinh </>, however 
the behavior with these potentials is very similar to the 
obtained with the potentials mentioned before. 

In Fig. [2] we show the accretion of the scalar field 
for all the potentials, in particular, we show the mass of 
the apparent horizon and the Misner-Sharp masses. The 
finding is that despite of the potential used, the final mass 
of the apparent horizon is nearly the same in all cases, 
that is, the potential has no effect on the final state of 
the resulting black hole. 

For both classes of configurations, positive and nega- 
tive ADM mass, we found that the apparent horizon mass 



decreases and converges asymptotically to the Misner- 
Sharp mass. As a particular example, we show in Fig. [2] 
the behavior of the apparent horizon mass, the Misner- 
Sharp mass and the even horizon location for the expo- 
nential potential case. 

In order to make sure that the decrease of the black 
hole is real, we fine-tuned a bundle of outgoing null rays 
for the positive and negative ADM mass cases, and esti- 
mated the location of the event horizon, since it is gauge 
invariant. We show the results in Fig. [5]where we present 
the location of the event horizon for the exponential po- 
tential case. The behavior shown in Fig. [5] is generic 
and applies to all the cases studied, that is, the event 
horizon -not only the apparent horizon- shrinks and is 
located approximately at the apparent horizon surfaces 
when the system has stabilized. 

Finally, in order to validate our numerical results we 
show in Fig. [3]the convergence of the L2 norm of the con- 
straints to zero as the resolution is increased with fourth 
order for two representative cases with the exponential 
potential. 

In Fig. Q] we summarize the results related to the ac- 
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crction mass rate and final mass of the apparent horizon 
for all various potential used. In each case we fine-tuned 
different scalar field profiles for the different potentials 
in order to use only one value of the ADM mass for the 
positive case and one value for the negative case. 

IV. CONCLUSIONS 

We present the full non-linear spherical accretion of a 
phantom scalar field with various potentials into a black 
hole. 

This accretion process is very efficient at reducing the 
black hole horizon area, which is a potentially important 
process that may impose restrictions on the existing black 
hole masses, or for instance primordial black holes. In 
our simulations we were able to reduce the mass of the 
initial black hole up to 50%. Attempts to reduce the area 
even further presents difficulties with our approach due 
to the steep gradients in the equations near the origin, 
and other techniques might be implemented in order to 
sort this problem out. 

We show that the final mass of the black hole's horizon 
is nearly independent of the potential used. This results 
is physically interesting because the process of area re- 
duction will apply independently of the cosmologically 
motivated potential used. 



We show that even though the Lagrangian (JfJ of the 
scalar field corresponds to a phantom field, at local scales 
the radial equation of state is always p/p > —I, and it 
depends on space and time, with values ranging from a 
cosmological constant equation of state lu = — I up to a 
stiff mater lu = I. 

We want to stress the importance of the equation of 
state of the scalar field given by the Lagrangian ((T|) , and 
point out that perhaps the big rip scenario would strongly 
depend on the homogeneity and isotropy of the scalar 
field. We have shown how much the equation of state 
can depend on space and time for the model of a scalar 
field, at least in the strong gravitational field regime. 

The astrophysical possibility of this accretion process 
strongly depends on the chance of having a phantom 
scalar field located near a black hole with a rather sharp 
profile. 
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FIG. 2. In the top panels we present the behavior of the different masses versus the coordinate time for the exponential 
potential and both, positive and negative ADM mass cases. The solid line represent the mass of the apparent horizon, the bold 
dashed line is the Misner-Sharp mass and the thin dashed line is the ADM mass of the space-time. In the positive ADM mass 
the apparent horizon of the black hole reduces up to 30 percent of its initial mass whereas in the negative case it is reduced up 
to 50 percent. In the bottom panel we show a set of outgoing null rays that help estimating approximately the location of the 
event horizon. For comparison we also show the location of the apparent horizon: at early stages of the accretion the apparent 
horizon lies outside of the event horizon, which is not surprising given that the energy conditions are not fulfilled and we notice 
that after the system has stabilized the event horizon and the apparent horizon coincide. 
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FIG. 3. We show the convergence of the Li norm of the violation of the Hamiltonian, momentum and G r constraints to zero 
for positive and negative ADM mass cases. 
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FIG. 4. In these plots, we show the behavior of the apparent horizon, Misner-Sharp and ADM masses with three self-interaction 
potentials: i) curves with subindex 1 correspond to the exponential potential, ii) curves with subindex 2 correspond to the 
quadratic potential, iii) Curves with subindex 3 correspond to the quartic potential. In the case of the exponential potential 
the decrease is faster than in the other cases, however the final mass of the horizon is not quite different considering that the 
scalar field profile is so free that the in-going profile depends on the scalar field potential. 



